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ABSTRACT 

An optical light curve of SUUMa type dwarf nova V1504Cyg taken by Kepler was 
analysed in order to study fast optical variability (flickering). We calculated power 
density spectra and rms-flux relations for two different stages of activity, i.e. quiescence 
and regular outbursts. A multicomponent power density spectrum with two break 
frequencies was found during both activity stages. The rms-flux relation is obvious only 
in the quiescent data. However, while the collection of all outburst data do not show 
this variability, every individual outburst does show it in the majority of cases keeping 
the rms value approximately in the same interval. Furthermore, the same analysis was 
performed for light curve subsamples taken from the beginning, middle and the end of 
the supercycle both for quiescence and regular outbursts. Every light curve subsample 
shows the same multicomponent power density spectrum. The stability of the break 
frequencies over the supercycle can be conhrmed for all frequencies except for the 
high break frequency during outburst, which shows variability, but with rather low 
conhdence. Finally, the low break frequency can be associated with the geometrically 
thin disc or its inner edge, while the high break frequency can originate from the inner 
geometrically thick hot disc. Furthermore, with our statistical method to simulate 
flickering light curves, we show that the outburst flickering light curve of V1504Cyg 
needs an additional constant flux level to explain the observed rms-flux behaviour. 
Therefore, during the outbursts another non-turbulent radiation source should be 
present. 

Key words: accretion, accretion discs - turbulence - stars: individual: V1504Cyg - 
novae, cataclysmic variables 


1 INTRODUCTION 

Dwarf novae (DNe) are a subclass of interacting binaries 
with a white dwarf as a primary star and a main sequence 
star as a secondary. The latter is Hlling its Roche lobe and 
mass is transfered to the white dwarf. With or without weak 
magnetic fields of the primary, a geometrically thin opti¬ 
cally thick accretion disc is formed. This class of interacting 
binaries is kn own as cataclysmic variables (CVs) (see e.g. 
IWarnerlll995l for a review). The accretion disc is the main 
source of variability on large intervals of time scales and 
energies. DNe exhibit typical alternating quiescent and ac- 
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tive (outburst) stages. The driving mechanism of this alter¬ 
nation is the viscous-thermal disc instab ility, generated by 
varia t ions in ionization state of hydroge n (lQsakJ|l974, iHoshil 
Il979l . [Meyer fc Mever-Hofmeistenll98ll . lLasotall2001f) . Dur¬ 
ing quiescence the mass accretion rate through the disc is 
much lower than during outbursts. While the mass is driven 
inwards, the angular momentum is transported outwards. 
Therefore, the different activity stages have different angu¬ 
lar momentum redistribution and disc structure. Different 
typical time scales of variability patterns are expected. 

An underlying accretion process is usually manifested 
in fast stochastic variability (sometimes also called flick¬ 
ering). Such oscillations are observed in optical and X- 
ray light curves of active galactic nuclei. X-ray bina- 
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ries and CVs. This variability reflects as red noise in 
power density spectra (PDS). The shape of such PDS can 
be comp licated by several characteristic break frequencies 
(see e.R. Sunvaev fc RevnivtsevI l200Cll . IScaringi et alJ l2012l . 
iDobrotka et alJ 20l4 l. The best example up to date in CVs 
is detected in high-cadence optical Kepler data of the nova- 
like system MV Lyr, t aken almost continuonsly over 600 days 
dScaringi et al.l[2012l l. The PDS show 4 different break fre¬ 
quencies. 

Second and even more fundamental characteristic of 
flickering is the correlation between variability amplitude 
and flux over a wide-range of timescales. This rms-flux re¬ 
lation has a typical linear trend simply meaning that the 
higher the flux, the lar ger the degree of variability. Standard 
shot-noise models (e.g. Terrelll [19721 1 with constantly decay¬ 
ing superposed flares (so-called stationary) seek to explain 
the broa dband noise variability incl uding the linear rms-flux 
relation. lUttlev fc McHardvl ll200ll ~) showed that the broad¬ 
band variability of the black hole binary Cyg X-1 and the 
accreting millisecond pulsar SAX J 1808.4-3658 are intrinsi¬ 
cally non-stationary, and as a consequence the rms scales 
linearly with the flux. Therefore, this linearity is a funda¬ 
mental characteristic of flickering activity. Following further 
studies of stochast ic X-ray variability of X-ray binaries and 
active galaxies bv lUttlev et al.l (l2005ll . the linear rms-flux 
relation implies that short and long variability time scales 
are coupled multiplicatively, and the flux of stationary data 
should have a log-normal distribution. 

This rms-flux linear relation is clearly a ubiquitous 
featu re in X-ray binaries and active galaxies jHeil et al.l 
[2 oT 2). But to detect such variability in CVs, long enough 
observations are needed because of longer variability time 
scales. This is now possible with the Kepler mission. So 
far, the rms-flux relation in CVs wa s detected in three 
systems, M VLvr llScaringi et al.ll2012ll . KIC8751494 and 
VI504 Cyg dVan de Sande et al.l 2015lb The most promising 
physical model describing the rms-flux relation are varia¬ 
tions in the accretion rate that are produced at different disc 
radii iLvubarskiil 19971. Kotov et al.l [20oi] . [Arevalo fc UttlevI 


I 2 OO 6 I . Ingram fc van der Klisll2013ll . This statement agrees 
with the recen t mod elling of the MV Lyr flickering by 
[Pobrotka et al.l (l2015l L The authors used their statistical 
method based on an u nstable turbulent mass accretion rate 
dPobrotka et al.ll2010h . and found that the linear rms-flux 
relation with log-normally distributed flux of stationary data 
is an inherent characteristic of their method. 

V1504Cyg is a member of a specific subset of the DNe 
showing also considerably larger superoutbursts that oc- 
cur in cycles , in addition to the regular outbursts (see e.g. 
IWarnerlll995l .lL aso t^l200ll for a review). The physical origins 
of these two types of outbursts are not the same. Two differ¬ 
ent physical scenarios have been proposed as explanations 
of the superoutbur st, i.e. enhanced mass transfer or tidal 
thermal instability llSchreiber et al.ll2004l l. The former sce¬ 
nario suggests that the superoutbursts are generated when 
the disc mass exceeds a critical value, while the latter is 
based on the outer disc radius expanding to a certain crit¬ 
ical radius (3:1 resonance radi us) where the t i dal a ctivity 
is trigerring the superoutburst. lOsaki fc Katd il2013[l anal¬ 
ysed the same data as we use in this paper, and studied 
the appearance of superhumps (a variability connected to 
superoutbursts). The authors concluded that the superout¬ 


burst was initiated by a tidal instability (as evidenced by 

the growing super hump). _ 

Finally. fVan de Sande et al.l (l2015l l studied the rms-flux 
relation of V1504Cyg during quiescence using the same 
data. The purpose of our work is to expand the analysis 
to the outburst stage and investigate the behaviour between 
two successive superoutbursts, i.e. during the supercycle. Af¬ 
ter the presentation of the analysed observations in Section[2] 
we present studies of the power density spectra (Section [3]l 
and the rms-flux relation (Section |4]). The results are dis¬ 
cussed and summarized in Sections [5] and [6] respectively. 


2 DATA 

The NASA Kepler mission offers a unique opportunity to 
study fast optical variability of a variety of accreting objects. 
The spacecraft was continuously monitoring a 116 square 
degree field-of-view with a cadence of 58.8 s until a second 
reaction wheel failure occurred. Such an almost continuou^ 
observation of the same sky region yields high quality optical 
light curves of several CVs over a time span of hundreds 
of dayfl This allows detailed studies of the power density 
spectra (PDS) over a wide frequency range. 

For our detailed light curve analysis we selected one of 
the longest light curves of a CV, i.e. the SU UMa DN system 
V1504Cyg (KIC number 7446357, orbital period of 1.67h). 
The light curve duration of about 1400 days consists of qui¬ 
escent emission with 118 fully covered regular outbursts dis¬ 
tributed into 11 full supercycles. The fast stochastic variabil¬ 
ity and the orbital pattern are obvious in both quiescence 
and outburst. 

We first extracted separate light curves from periods 
of quiescence and outburst. We first identified all outburst 
peaks with a simple algorithm. Subsequently, we determined 
the times when the rise to outbursts ended and when the 
decline from the outburst started by taking the time deriva¬ 
tives of flux and applying certain limits to the derivatives af¬ 
ter and before the peak (derivative limit = 1000 electrons/s). 
Outburst light curves were extracted from the times in be¬ 
tween these time intervals. Time intervals of quiescence be¬ 
tween outbursts were identified by applying an upper flux 
limit (dashed curve in Fig.[T} in between the times of decline 
of one outburst and rise into the respective next outburst 
again determined by a derivative limit (derivative limit = 
1500 electrons/s). The limits on the derivatives and flux were 
chosen visually. Both light curves with details are shown in 
Fig. [ 1 ] and all starting and ending times of all individual 
light curve subsamples are summarized in Table [T] 

Subsequently we selected the start, middle and end 
stage of the supercycle in order to study any evolution for 
both quiescence and outbursts. Therefore, we collected all 
quiescent intervals lasting as first (later called start), sec¬ 
ond (start-fl) and third (start-|-2) after the superoutburst. 
A similar selection was done with quiescent intervals before 
the superoutburst, i.e. we collected last (end), second (end- 
1) and third (end-2) interval before the superoutburst. As 


^ The satellite performed quarterly 90 deg rolls in order to have 

the solar panels face to the Sun. _ 

^ For Kepler CVs light curves summary see IVan de Sande et al.l 
ll2015h 
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Figure 1. Analysed light curves. Upper left panel: all data divided into quiescent (below the dashed line) and outburst subsamples 
(above the dashed line). Lower left panel: supercycle example with marked supercycle subsamples, i.e. 1 - start, 2 - start+1, 3 - start+2, 
4 - middle, 5 - end-2, 6 - end-1 and 7 - end (see text for details). Right panels: examples of quiescent and outburst intervals. The dots 
show the original light curve and circles are the selected data (outburst or quiescence). 


middle of the supercycle we choose the closest interval to the 
half time of the supercycle. The same was performed with 
the outburst intervals (inset of Fig. [TJ. 


3 PDS ANALYSIS 

For the PDS a nalysis we applied the Lomb-Scargle algorithm 
llScargld Il982l ) which can handle gaps in the light curves. 
The low frequency end of the studied PDSs usually depends 
on the light curve duration. In order to study differences 
between outburst and quiescence both PDSs must be com¬ 
parable on the studied frequency interval. Therefore, we set 
this limit empirically (see next Section). The high frequency 
end is usually limited by the white noise or power rising of 
the PDS to the Nyquist frequency. We empirically defined 
this limit to Hz. The quiescent data were divided into 

one day intervals. Every quiescent interval yields an individ¬ 
ual PDS. The outburst data were not divided into shorter 
subsamples. Every outburst light curve was used for indi¬ 
vidual PDS calculation. Finally, a mean PDS was calculated 
from all individual PDSs. We subsequently binned the mean 
PDS into equally spaced bins and the mean values with the 
errors of the mean was fitted with a broken power law. 

3.1 All data 

Fig. El shows (left column) the PDSs of all quiescent and 
outburst data. The orbital frequency and its first harmonic 
can clearly be seen and are marked by dashed lines. While 
the quiescent PDS shows a significant peak around / = 
Hz (see the low resolution insets in Fig. [2]), the out¬ 
burst PDS is influenced by the overall inherent outburst 
shape in this low frequency region and shows a steep rise 
toward lowest frequencies. A comparable behaviour of both 
PDSs starts above the orbital frequency. Therefore, in order 
to focus our study on the common behaviour of both PDSs, 
we choose the frequency / = 10“®'® Hz as a low frequency 
end of the studied PDSs. 


The multicomponent shape of the PDS with two break 
frequencies /l (low frequency) and /h (high frequency) is 
clearest in the quiescence PDS. The outburst PDS has one 
clear break frequency while a multicomponent shape with a 
second break frequency around / = 10“^'^ Hz may be pos¬ 
sible. We fitted the binned PDSs with broken power law 
modelu. In the case of outburst we used a multicomponent 
modeO, while in the quiescent PDS two separate broken 
power law fits were more adequate due to the convex shape 
between the two broken power laws. The resulting fitted 
PDS parameters with 1 -ct errors (derived from test) are 
summarized in Table[2] In the multicomponent fit of the out¬ 
burst PDS, we had problems with the fitting procedure while 
keeping /h constant in order to get the curve. The fit 
often converged to deformed unacceptable solutions. There¬ 
fore, in order to derive the l-cr error of the outburst /h, we 
fitted a single power law model to the selected PDS interval. 
Finally, the x^ residuals of the fits often yield asymmetric 
curves. Therefore, in TableEJwe show lower and upper limits 
for clarity, instead of plus/minus errors. 

In general, we can conclude that the multicomponent 
shape of the PDS with two characteristic break frequencies 
is present in both activity stages. Furthermore, /l is sig¬ 
nificantly different, while /h agree within the l-cr errors for 
both activity stages. 


® The bin containing the orbital harmonic frequency was ex¬ 
cluded from the fitting procedure. 

We reduced also the fitting frequency interval because of large 
scatter in the data in the lowest frequency region which influenced 
the iteration procedure in some cases. We used / = 10“®'^ Hz as 
the low-frequency end instead of / = 10“®'® Hz. This does not 
influence the study, because the omitted low-frequency region is 
useless and the searched break frequency is clearly around / = 
10-3 ° Hz. 










4 A. Dobrotka 


Table 1. Starting and ending dates of individual light curve intervals for all selected quiescence and outburst subsamples. The number 
of full supercycles is 11 which yields 11 middle light curve intervals, while part of the supercycle at the beginning and at the end of the 
light curve yields 12 light curve intervals in the case of other subsamples. Where an observing gap is interrupting the light curve, or the 
searched interval is somehow not distinguishable, the light curve interval is missing. 


supercycle 

quiesc. 

quiesc. 

outb. 

outb. 

supercycle 

quiesc. 

quiesc. 

outb. 

outb. 

position 

start 

end 

start 

end 

position 

start 

end 

start 

end 


(days) 

(days) 

(days) 

(days) 


(days) 

(days) 

(days) 

(days) 


start 256.3280 

381.4602 
491.8608 

586.3605 
696.4608 

817.6604 

930.1603 
1057.0605 

1173.1603 
1295.6785 

1534.9608 

start+1 262.6602 

388.3606 
499.0601 
594.1606 
704.8580 
822.7608 
936.9335 
1064.2610 

1180.3606 
1302.4604 
1419.1587 
1542.7584 

start+2 273.7607 

506.8611 

607.0575 

716.8602 

828.1603 
945.1608 
1072.3456 
1188.4607 
1311.4603 
1429.9607 
1556.2606 

middle 306.4605 

424.9608 
528.7606 

627.1603 
736.0604 

862.6604 


258.4775 

260.8069 

385.6599 

386.3908 

496.3605 

496.9292 

591.7604 

592.2522 

702.1602 

702.7916 

820.3605 

821.3372 

934.6599 

935.2852 

1061.5598 

1062.2777 

1177.9617 

1178.4786 

1300.0602 

1300.6895 

1540.0604 

1540.7157 

271.0581 

271.7412 

393.4601 

394.1099 

504.6044 

505.0179 

604.0599 

604.8896 

713.8601 

714.7352 

825.7599 

826.5732 

942.4601 

943.2925 

1070.0155 

- 

1185.7600 

1186.5957 

1308.7603 

1309.5906 

1426.9598 

1427.7765 

1553.5598 

1554.1177 

280.5307 

281.6477 

- 

402.5966 

514.3598 

514.9804 

623.5601 

- 

719.5580 

- 

831.1600 

832.2969 

949.3600 

950.2032 

1077.7603 

1078.4399 

1195.9601 

1196.7468 

1318.3601 

1319.0753 

1436.3187 

1437.0366 

1560.7616 

1561.5885 

313.6603 

314.5675 

431.8600 

432.8933 

535.9603 

526.3112 

646.3605 

624.4659 

753.7604 

754.5893 

868.0599 

861.0216 


261.1876 end-2 

386.7041 

497.2998 

592.6248 

703.1669 

821.7221 

935.6638 

1062.6557 

1178.8008 

1300.9974 

1541.2333 

271.9762 end-1 

394.5485 

505.3700 

605.2485 

714.9627 

826.8824 

943.6842 

1186.9240 

1309.8256 

1428.1694 

1554.5019 

281.9808 end 

402.9849 

515.3236 


832.7492 

950.6010 

1078.8309 

1197.1276 

1319.4240 

1437.4255 

1561.9712 

314.9619 middle 

433.2836 

526.6729 

624.8167 

754.9625 

861.3267 


215.8605 

220.9621 

339.7604 

346.9603 

444.4605 

454.3602 

782.5609 

786.1599 

888.4603 

898.3601 

1030.3605 

1033.3600 

1264.3604 

1267.9602 

1501.9603 

1506.1602 

224.5604 

233.2603 

352.4397 

359.5599 

458.2604 

464.5603 

539.5853 

553.9600 

649.9608 

667.9600 

790.0607 

793.3602 

902.2603 

905.9364 

1037.2606 

1040.5604 

1143.1606 

1147.0601 

1270.3605 

1272.4604 

1390.9608 

1394.2600 

1509.7605 

1513.3602 

236.8606 

241.6584 

363.1602 

368.2603 

468.4606 

478.6601 

558.1605 

573.4600 

672.1602 

683.5603 

796.3605 

800.26001 

909.4606 

916.9599 

1044.4604 

1045.0604 

1148.8602 

1152.7556 

1276.0603 

1281.4603 

1396.3605 

1398.1586 

1516.9606 

1522.0602 

979.3605 

987.7599 

1105.9603 

1110.1566 

1219.3602 

1224.4600 

1342.6604 

1349.2604 

1467.1606 

1472.5604 


213.5107 

213.8853 

337.1763 

337.557 

442.0567 

442.2127 

779.9269 

780.3152 

886.2139 

886.7254 

1027.9310 

1028.2879 

1133.0365 

1133.4343 

1262.0771 

1262.8127 

1371.3064 

1371.3323 

1499.0002 

1499.3572 

222.3209 

222.6009 

348.2367 

348.6249 

455.2695 

455.7606 

536.9200 

537.2878 

647.3603 

647.7540 

787.5201 

788.1236 

899.4737 

899.9335 

1035.0144 

1035.6015 

1140.9021 

1141.6622 

1268.8763 

1269.1440 

1388.1697 

1388.5763 

1507.2479 

1507.7710 

234.4168 

234.8459 

360.5455 

360.9752 

466.1433 

466.5820 

555.1349 

555.4414 

669.0184 

669.3691 

794.0964 

794.4377 

907.2884 

907.4853 

1041.7701 

1042.5397 

1157.7912 

1158.3014 

1274.1299 

1274.2362 

1394.8117 

1395.0637 

1514.3703 

1515.0671 

988.5288 

988.9423 

1110.9556 

1111.6421 

1225.6397 

1226.0511 

1350.0273 

1350.4074 

1465.0343 

1465.4021 


3.2 Supercycle subsamples 

We performed the same analysis for other data subsamples 
collected for different supercycle phases (11 full covered su¬ 
percycles). Fig. |3] shows low and high frequency parts with 
the broken power law fits separately for quiescence, and 
Fig. [4] shows the outburst PDSs with the multicomponent 
fits. Clearly /l and /h are present in all supercycle phases 
in both cases. The outburst PDSs show considerably larger 
scatter than the quiescence data, likely because of much less 


data. This is visible mainly in the last (end) PDS where it 
was difficult to find a best fit. The result is not robust against 
some factors, i.e. initial parameter estimate or the frequency 
interval over which the fit was performed. We show two so¬ 
lutions, i.e. the multicomponent fit as it converged, and vi¬ 
sually motivated simple broken power law. Furthermore, in 
the case of outburst /h the residuals calculations did not 
yield the searched l-cr confidence limits in every case. Con¬ 
siderable reductions of the confidence limit are needed to 
determine the error from the residuals curves but such 
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Figure 2. Left column: mean PDS of quiescent and outburst data. Solid lines are the broken power law fits. The dashed lines shows 
the orbital frequency / = Hz (1.67 h) with its first harmonic / = Hz (0.84 h). The insets show low resolution mean PDSs 

with larger x-axis range. Right column; Rms-flux relations calculated from both quiescent and outburst data. The dashed lines are linear 
fits. 


Table 2. Fitted low and high (/l and /h) break frequencies and a coefficient from the rms-flux linear fit for different subsamples with 
l-(7 limits. ’’Stage” is the activity stage following the disc instability model and ’’sample” is the analysed data subsample. Two solutions 
are listed for the subsample ’’outburst, end” because it is model dependent. The outburst l-tr limits within parentheses are derived 
from single broken power law fits over a reduced frequency interval, instead of the original multicomponent fit performed over the full 
frequency interval. 


stage 

sample 


log(/L) 



log(/H) 


a 

quiescence 

all 

-3.37 

-3.38 

-3.39 

-2.39 

-2.41 

-2.42 

0.030 ± 0.002 

outburst 

all 

-2.94 

-3.00 

-3.03 

(-2.27) 

-2.33 

(-2.40) 

0.0017 ± 0.0003 

quiescence 

start 

-3.36 

-3.39 

-3.42 

-2.35 

-2.39 

-2.45 

0.030 ± 0.003 

quiescence 

start+1 

-3.36 

-3.41 

-3.46 

-2.38 

-2.40 

-2.43 

0.029 ± 0.003 

quiescence 

start+2 

-3.34 

-3.42 

-3.46 

-2.40 

-2.42 

-2.44 

0.046 ± 0.004 

quiescence 

middle 

-3.26 

-3.30 

-3.42 

-2.34 

-2.39 

-2.43 

0.056 ± 0.003 

quiescence 

end-2 

-3.32 

-3.37 

-3.40 

-2.35 

-2.39 

-2.43 

0.029 ± 0.004 

quiescence 

end-1 

-3.38 

-3.40 

-3.45 

-2.33 

-2.42 

-2.45 

0.044 ± 0.003 

quiescence 

end 

-3.33 

-3.43 

-3.47 

-2.37 

-2.40 

-2.44 

0.050 ± 0.008 

outburst 

start 

-2.97 

-3.07 

-3.18 

(-2.39) 

-2.43 

(-2.44) 

0.001 ± 0.002 

outburst 

start-hl 

-2.83 

-2.94 

-3.16 

- 

-2.33 

- 

0.002 ± 0.001 

outburst 

start-|-2 

-2.79 

-2.92 

-3.03 

(-2.44) 

-2.47 

(-2.51) 

0.003 ± 0.001 

outburst 

middle 

-2.72 

-2.90 

-3.10 

- 

-2.52 

- 

0.0018 ± 0.0006 

outburst 

end-2 

-2.73 

-2.95 

-3.10 

(-2.28) 

-2.32 

(-2.34) 

0.0017 ± 0.0005 

outburst 

end-1 

-2.72 

-2.82 

-3.14 

(-2.33) 

-2.41 

(-2.44) 

0.0022 ± 0.0004 

outburst 

end 

- 

-2.59 

- 

- 

-2.29 

(-2.46) 

0.0018 ± 0.0005 



- 

-3.01 

- 






low confidence intervals are meaningless. Therefore, some 
fitted /h values are questionable, and we can only conclude 
the presence of /h but with a highly uncertain value. All 
fitted break frequencies with derived 1-cr limits are again 
summarized in Table [2] 


Fig.[5]shows the evolution of the fitted break frequencies 
as a function of the supercycle phase. All frequencies except 
the outburst /h are stable within the derived l-cr confidence 
limits, and no trend can be deduced. Just in the case of out¬ 
burst /l a rising trend is possible but due to large errors 
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Figure 3. 

shows the 


Power density spectra of quiescent data with broken power law fits for different supercycle stages. The vertical dashed line 
first harmonic of the orbital frequency. 


still speculative. The outburst /h values are scattered, and 
some values do not agree within the l-cr confidence inter¬ 
val. This suggests non-stability during the supercycle. But 
slightly larger confidence limits than the used l-cr would 
yield a match within the confidence frequency limits. There¬ 
fore, the instability of the outburst /h is possible with rather 
low probability of approximately l-cr. 


3.3 X-ray data 

V1504Cyg was in the field of an XMM-Newton observation 
taken 2014 October 15 00:45-07:08UT (ObsID 0743460201) 
that, at the time of writing is still proprietary. With the kind 
permission of the PI, J.H.M.M. Schmitt, we extracted the 
X-ray light curve of V1504 Cyg from the two MOS detectors 
(the source was not within the field of view of the pn) in 50 s 


binnin g and studied it in the same way as in iDobrotka et ahl 
(l2014h but only found white noise. This may be due to the 
count rate of < 0.1 counts per second in each MOS light 
curve, 100 times lower than in RU Peg which we tested by 
simulations. We took a subsample of the Kepler data and 
reduced the count rate to the same level as observed with 
XMM-Newton. Subsequently we added Poisson noise to the 
simulated light curve in order to get a comparable flux his¬ 
togram than the X-ray data. The break frequencies are still 
in the simulated light curve, but can not be recovered, sup¬ 
porting our interpretation that the count rate is too low. 
Thus, if the break frequencies are also present in the X-ray 
light curve, they would be buried in the noise. 

While re-binning can increase the signal to noise ratio, it 
did not change the negative results. Furthermore, the larger 
bin size reduces the time resolution of the light curve and any 
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Figure 4. Same as Fig. [3] derived from the outburst light curves. 


binning larger han 100 s is already larger than the required 
time resolution for a detection of /h. We thus conclude that 
V1504 Cyg is too faint in X-rays to complement our optical 
study with an X-ray PDS with the available sensitivity of 
the current instrumentation. 


4 RMS-FLUX RELATION ANALYSIS 


4.1 Data 

The absolute rms amplitude of variability is defined as 
square-root of the variance, i.e. 




N - 


iV 


( 1 ) 


where N is the number of data points, Xi is the i-th point 
and X is the mean value of all Xi. We subdivided all light 
curves into small parts each containing 10 points (we tried 
more values, i.e. 25, 50, 100 and 200) of the light curve. 
Subsequently for each small part we calculated the corre¬ 
sponding rms and mean flux. The binned data by flux were 
fitted with a linear function. 

The right column of the Fig. [2] shows the rms-flux 
relations of all quiescent and outburst data. The qui- 
escent rms-flux ca n be clearly seen as expected from 
IVan de Sande et al.l (120151 '). but the outburst case does not 
show the same characteristics. The linear fit shows a rising 
trend, but over a wide flux range. The parameter a of the 
linear function rms = aijj + b {a is the gradient, b is the verti¬ 
cal offset and i/i is the flux) is only 0.002 compared to 0.03 in 
quiescence. Furthermore, the rms-flux data are always more 
scattered toward the highest flux, which can generate an 
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Figure 5. Evolution of the fitted PDSs and rms-fiux parameters during the supercycle. The error bars are l-cr errors. When this 
confidence limit was not reached in the test (confidence interval lower than the selected l-c), a vertical line over the whole y-axis 
interval is used instead. 


artificial linear rising trend. We selected the less scattered 
data region with flux between 15000 and 20000 electrons/s, 
which yields a = —0.001. Clearly it is hard to talk about a 
rising linear trend which should be robust against flux range 
selection as is obvious in the quiescent case. 

We performed the same procedure for all other data 
subsamples and the resulting rms-flux relations are shown 
in Fig. [6] The basic difference between quiescent and out¬ 
burst rms-flux is robust against the supercycle phase. The 
evolution of a is shown in Fig. [S] (rms calculated from 10 
and 50 points is shown) and all a values are listed in the last 
column of table [2] The quiescent a values are scattered in 
the supercycle without any significant trend, and the ’’evo¬ 
lution” shape is not robust against the selected number of 
data points over which the rms is calculated. All outburst 


a values are still of the order of 0.001 (0.003 as maximum). 
Therefore, as for all outburst data, any linear rms-flux trend 
is not convincing. 


In the last column of the Fig. |6]we show linear fits to 
the rms-flux data of individual outbursts. Clearly, while the 
typical rms-flux relation is hardly detectable in the collection 
of all outburst data, the individual outbursts do show this 
relation in most cases. Moreover, linear fits of all individual 
outbursts lie approximately between 50 and 150 electrons/s 
(the rms-flux values are concentrated between 50 and 200 
electrons/s). 
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Figure 6. Rms-flux relations with corresponding linear fits for quiescent and outburst data. The right column shows linear fits to 
individual outbursts. 


4.2 Simulations 

In lPobrotka et all (|201Gh we developed a statistical method 
to simulate light curves dominated by red noise, also known 
as flickering. The method is based on the simple idea of dis¬ 
crete angular momentum transport between two adjacent 
concentric rings in a geometrically thin, optically thick, hot 
ionized steady-state disc. The angular momentum is car¬ 
ried by discrete turbulent bodies with an exponential dis¬ 
tribution of the dimension scales. When a turbulent body 
penetrates from one ring to the adjacent one, it changes 
the tangential velocity and distance from the center, while 
the mass is conserved. This changes the angular momen¬ 
tum. The sum of changes in angular momentum of all bodies 
must be equal to the total difference in angular momentum 
between two concentric rings. The distribution function of 
dimension scales is equal to the distribution function of tur¬ 


bulent event time scales because the dimension scale divided 
by the local viscous velocity gives the time scale of a tur¬ 
bulent event. All events liberate energy in the form of a 
flickering flare. The distribution functions between pairs of 
rings are summed over the entire disc (from the inner disc 
radius to the outer disc radius), which yields the final dis¬ 
tribution function of flare time scales. Following this distri¬ 
bution function, flickering flares are randomly redistributed 
into a synthetic light curve of the same duration and sam¬ 
pling as the observations. Such synthetic light curves are 
subsequently analysed. 

Furthermore, in iDobrotka et al.l l2015l we applied the 
simulation method to the flickering study of MVLyr sys¬ 
tem observed by Kepler. We showed that the simulated light 
curves show the typical linear rms-flux relation. MV Lyr is a 
nova-like system with the disc in the so-called hot stage with 
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fully ionized hydrogen fsee lLaso^ ll200ll ') for a review). This 
hot stage, together with the steady state, is the essential as¬ 
sumption of our model. The studied data of V1504Cyg in 
this paper do not satisfy these basic assumptions, i.e. the 
quiescence is a cold stage with no ionized hydrogen, and the 
outburst is not in a steady state. Therefore, we do not use 
these simulations to study the PDS in detail, but we use it 
to study how it influences the rms-fiux relation. The latter 
is based on flickering flare redistribution and the principle is 
robust against the underlying physical model of the accre¬ 
tion, i.e. non steady/steady state, quiescence/out burst etc... 

It appears that using a larger number o f flares per light 
curve increases the rms and also the flux ( Dobrotka_et_^ 
2015) while the PDS remains unchanged IIDobrotka et al. 
20141 . In order to investigate how this plays a role in dif¬ 
ferent rms-fiux behaviour in quiescence and outburst of 
V1504 Cyg, we simulated different light curves with red noise 
(Fig.0. As input parameters we used the known orbital pe¬ 
riod and typical values in CVs for the rest (1 Mq as primary 
mass, 10^^ gs~^ as mass accretion rate through the disc, in¬ 
ner disc radius of 10® cm and outer disc radius estimated as 
0.9 times the Roche lobe radius jf). 

In panel A of Fig. [7] we show a simple simulated light 
curve with constant mass accretion rate and corresponding 
rms-fiux relation. This simulated light curve is equivalent 
to the observed quiescent phase of V1504Cyg shown in the 
upper panel of Fig. [2] The typical rms-fiux relation is clear 
and corresponds to the observed relation. 

In panel B we constructed a rising light curve by adding 
small light curve subsamples one by one with increasing 
mass accretion rate and proportionally higher numbers of 
flarefl The basic idea is based on the assumption that the 
mass accretion rate scales linearly with the transported an¬ 
gular momentum and thus the number of flickering flares. 
We gradually changed the mass accretion rate from 1 to 
9 X 10^^ gs“^ and the number of flares from 10000 to 90000. 
Clearly, the resulting light curve shows the typical linear 
rms-fiux relation. This si mulation is equivalent to the rising 
flux of MVLyr studied bv IScaringi et all (l2012fl . 

Panel C is the same as panel B, without increasing the 
number of flares. The rms-fiux relation is different from the 
two previous cases and is more similar to the observed out¬ 
burst case of V1504 Cyg shown in lower panel of Fig. [5] 

In panel D we simulate the outburst light curve with 
typical gaps between the maxima. We made the same sim¬ 
ulations as in panel B, but we redistributed the light curve 
subsamples randomly to mimic the behaviour of the ob¬ 
served data in Fig. [T] The resulting rms-fiux relation is 
present as expected. 

The light curve in panel E is generated in the same way 
as in panel C, i.e. every light curve subsample is generated 
with the same number of flares. To get higher flux during 
the outbursts we added a randomly selected constant level 
to every small light curve segment. The resulting rms-fiux 
relation is similar to case C, and to the observed outbursts 
of V1504 Cyg. 


® Different parameter values do not influence the behaviour we 
want to study/show. The principle studied in this paper is based 
on light curve construction, not on real flickering statistics. 

® In our algorithm we can not change the parameters continually. 


For both simulated outburst light curves (panels D and 
E) we added a panel with linear fits to every individual out¬ 
burst light curve segment. Not every one, but the majority 
of the fits have the rising trend, which agrees with the obser¬ 
vations. The presence of decreasing linear fits in simulations 
or observations are due to short analysed data sets. In long 
enough light curves we never noticed such deviation (panel 
A). 


5 DISCUSSION 

In this paper we analysed Kepler d&ta of the DN V1504 Cyg, 
focusing on PDS and rms-fiux relation differences between 
quiescence and regular outbursts. Further, we investigated 
the evolution during the supercycle. 

5.1 PDS 

Both the quiescent and outburst PDS have a multicompo¬ 
nent shape with two break frequencies, with a lower and a 
higher frequency break, /l and /h respectively. It is believed 
that these break frequencies are the imprint of the under¬ 
lying physical processes an d disc structures llScaringi et al.l 
12012 . IPobrotka et al.l[2015h . We found discrepant values of 
/l for quiescence and outburst episodes, while /h agrees 
within the errors in both PDS. Therefore, /l should rep¬ 
resent a structure or process which is changing during the 
viscous-thermal instability while /h should arise from stable 

structures, unaffected by the outbursts. _ 

Following the disc instability model (see iLaso t^ I 2 OOII 
for a review) which describes the viscous-thermal instabil¬ 
ity, the mass accretion rate during the quiescence is too 
low causing inefficient cooling in the innermost parts of 
the disc. This is the reason for gas evaporation, and an 
optically thin geometrically thic k corona forms (see e.g. 
iMever fc Mever-HofmeistCTll994h . On the other hand, dur¬ 
ing the outburst and hot ionized steady state of nova like sys¬ 
tems, the mass accretion rate is strong enough, the density in 
the central parts ensures efficient cooling and the disc is (al¬ 
most) not truncated. One would expect that the hot X-ray 
corona should be abs ent, but a sandwiched model is a possi¬ 
bility as proposed bv IScaringi et all ll2012h . i.e. the geomet¬ 
rically thin disc is fully developed down to the white dwarf 
surface with a corona above the disc material developed from 
the white dwarf surface to a certain radius lower than the 
outer edge of the geometrically thin disc. Therefore, the in¬ 
ner radius of the geometrically thin disc should be different 
between quiescence and outburst, while the hot corona can 
be present during both stages. This suggests that the fre¬ 
quency /l is generated in the thin disc or its inner edge, 
while the frequency /h in the corona. A similar interpre¬ 
tation was proposed for RUPeg in quiescence where X-ray 
observations directly sug gest that the low freq uency is gen¬ 
erated in the inner disc (iDobrotka et al.ll20l4 fn. Moreover, 

^ We do not investigate the V1504Cyg PDS by simulations in 
this paper, because our method is not yet adapted to quiescent 
discs (Dobrotka, in preparation), nor to non-steady states during 
outbursts. Future work is planned to perform all three simulations 
together, i.e. flickering generated by 1) the whole quiescent disc, 
2) unstable mass accretion rate from the inner quiescent disc as in 
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Figure 7. Simulated light curves (left column) and corresponding rms-flux binned data. Middle column shows the rms-fiux data taken 
from full light curves with linear fits as dashed lines with a being linear coefficients of equation rms = aijj + b (where p is the flux). 
The two lower right panels show linear fits to individual simulated light curve subsamples (outbursts). The upper right square panel is 
a detail of the upper middle panel. Flux is in arbitrary units (au). A - constant mass accretion rate, B - increasing mass accretion rate 
with correspondingly increasing number of flickering flares per light curve subsegment, C - the same as B but with constant flickering 
flare number per light curve subsegment and an additional vertical offset is added instead, D and E - the same as B and C respectively, 
but light curve subsegments are randomly redistributed in order to simulate the outbursts. See text for details. 


the inner disc is moving inwards during the outburst where 
the typical time scale s are shorter (Dobrotka, in preparation, 
IPobrotka et al.l|2012l ). Therefore, a higher break frequency 
during outburst is plausible. 

Detailed investigations of the two break frequencies in 
quiescence did not show any significant variability during the 
supercycle. We can even say, that /h is stable with high con¬ 
fidence as the measurement error decreases with the investi¬ 
gated frequency. Therefore, the disc structure/shape should 
be stable during the supercycle, or hardly detectable in a 
PDS measurement. 

The situation in outburst is slightly different. While all 
values of /l agree within the errors during the supercycle, 
a rising trend is still possible. But the measurement error 
does not allow stronger conclusions. Furthermore, the fre¬ 
quencies /h do not agree within the calculated 1-a errors, 
which suggests a rather unstable coronal structure. This con¬ 
clusion is not surprising because we expect a corona to be 
present for low accretion rate d uring quiescence (see e.g. 
iMever fc Mever-Hofmeisterj|l99J ). while a high mass accre- 

RU Peg in iDobrotka et ^ ll2014h and 3) full disc during outburst 
(with the assumption that the disc is fully developed down to the 
white dwarf surface during the outburst). 


tion rate during the outburst suggests the opposite. But the 
corona is still possible in low density boundary regions above 
the disc photosphere like in the Sun, which can explain the 
sandwiched model. Finally, in order to strengthen our inter¬ 
pretation, direct X-ray detection of /h would be needed. As 
described in Section E3I the X-ray count rate in available 
XMM-Newton data is too low for meaningful conclusions. 


5.2 Rms-flux relation 

In iDobrotka et ^ (l2015l) we investigated the rms-flux rela¬ 
tion of the nova like CV MVLyr. We found that the rms- 
fiux relation is generated by a varying number of superposed 
flickering flares over the simulated light curve, because the 
the rms and the flux increase with the number of flares. Such 
local increases of the number of flares can be due to a random 
redistribution of the flares or by the larger mass accretion 
rate. The former is shown in simulations shown in panel A 
in Fig. El while the latter in panel B. The model used is not 
adequate for the purpose of this paper, but the principle 
of light curve construction is independent of the underlying 
physics. Panel A can be a representation of a quiescent light 
curve as shown in Fig.[T]with the corresponding rms-fiux re¬ 
lation shown in Fig.[2| Panel B has a rising trend with clearly 
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noticeable rising variability amplitude. This simulated light 
c urve resembles v ery well the Kepler data of MV Lyr in Fig. 1 
of IScaringil ll20l4l for intervals 15 0-450, 800-900 or 1000-1100 
days. While l^caringi et al.l il2012l l show the rms-flux relation 
only for a small light curve subsample, we show the relation 
for the full simulated light curve. But it is valid also for small 
subsamples. 

However, the MV Lyr system is a nova like CV with 
the disc in an ionized hot state. The same is valid for the 
outburst. Therefore, we expected a similar behaviour. Panel 
D of Fig. [7] shows an outburst light curve generated with 
this assumption. The discrepancy with observations is clear. 
In order to get the observed rms-flux relation in simulated 
outburst data, the light curve should be constructed in a 
different way. Panels C and E show the case with a constant 
number of flares but with added constant flux vertical offset. 
Therefore, the original assumption that is valid for the nova¬ 
like system MV Lyr (the mass accretion rate increases with 
the number of turbulent events) is needed to transport the 
enhanced angular momentum, is violated. Apparently both 
hot ionized states are not identical. Whether the disc is in 
a steady or non-steady state is essential while the activity 
stage does not matter. The physical reason for this differ¬ 
ent behaviour is not known to us, but at least it suggests 
that the quantity of flickering flares does not depend on the 
mass accretion rate during non steady-state outbursts while 
it does in steady-state discs of nova like systems. A non- 
turbulent steady radiation source during the outburst but 
not present or not dominant in the steady state of nova like 
CVs would resolve the puzzle. 

The disc instability model as a driving mecha nism of 
DN outbursts is described in detail in lLaso M ll200lll . Their 
Fig. 10 shows the surface density evolution during the out¬ 
burst process. Apparently the density in the innermost 
parts of the disc as a dominant flickering source (see e.g. 
iDobrotka et aklboid . iDobrotka et al.|[2012li increases by two 
orders of magnitudes during a relatively short time. This 
abrupt change in flow characteristics is probably the cause 
for significantly different physical conditions compared to 
the steady state accretion of nova like systems. Therefore, 
the radial dependence of the mass accretion rate in the disc 
is probably not the only difference between outbursting and 
nova-like discs, but also the hydrodynamical conditions are 
significantly different. The turbulent characteristics of the 
outbursting flow is somehow saturated yielding a very sim¬ 
ilar almost equal rms value irrespective of the total flux. 
Therefore, the study of the enigmatic viscosity in accretion 
discs should be oriented in three different regimes instead 
of two (cold not ionized, hot ionized stages), i.e. quiescence, 
outburst and steady state accretion of hot ionized disc in 
nova like systems. 

Finally the individual outburst light curves show the 
rms-flux relation in the majority of cases in both the simu¬ 
lations and observations. Some deviations from the typical 
rising linear relation are clear in both cases due to the short 
light curve duration. This rms-flux relation is generated by 
the random flare redistribution which is in agreement with 
the original assumption of underlying turbulent accretion. 


6 SUMMARY 

We analysed Kepler data of the dwarf nova V1504 Cyg. We 
searched for differences between optical flickering in quies¬ 
cence and outburst. The results can be summarized as fol¬ 
lows: 

(i) Both power density spectra show a multicomponent 
behaviour with two break frequencies. 

(ii) The low break frequency is higher in outburst data 
while the higher break frequency agrees within the l-cr er¬ 
rors. 

(iii) The break frequencies are stable during supercycle 
except for the high frequency in outburst light curves which 
show no significant trend. 

(iv) We conclude that the low frequency break origi¬ 
nates from regions affected by the outbursts while the sta¬ 
bility of the higher break frequency indicates an origin that 
is not affected by outbursts. 

(v) We propose that the low frequency break can be 
associated with the accretion disc or its inner edge, while the 
high frequency break can be generated within the inner hot 
geometrically thick corona. We are planning to study this 
model by detailed simulations of the power density spectra 
in a future work. 

(vi) The rms-flux relation is confirmed in quiescence 
but not in outburst when analysing all data as a single data 
set. Individual outburst light curves show the typical linear 
rms-flux relation. 

(vii) The same as previous point is valid for individual 
supercycle subsamples and no significant monotone evolu¬ 
tion trend of the flux gradient is found. 

(viii) We simulated both the observed quiescent and 
outburst rms-flux relations to test various scenarios. In the 
quiescent case a random redistribution of flickering flares 
in the simulated light curve is generating the rms-flux re¬ 
lation while for outbursts the basic idea was based on the 
assumption that the number of superposed flickering flares 
is proportional to the mass accretion rate. This we found 
valid for the hot disc in the nova-like system MV Lyr but 
not for hot discs in outburst. For the latter, the number of 
flickering flares must be kept unchanged, but a constant flux 
to the light curve must be added. 

(viii) While nova-like steady state discs and non-steady 
state discs in dwarf novae in outburst are both in a hot ion¬ 
ized state, the radiation generation mechanism is not iden¬ 
tical. The previous point suggests that flickering in both 
cases is generated by the turbulent mass accretion in a disc, 
but during the outbursts another additional significant non- 
turbulent source of radiation must be present. 
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